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Abstract 

The paper is devoted to the prospects of using the laser radiation interaction with plasmas in 
the laboratory relativistic astrophysics context. We discuss the dimensionless parameters charac- 
terizing the processes in the laser and astrophysical plasmas and emphisize a similarity between 
the laser and astrophisical plasmas in the ultrarelativistic energy limit. In particular, we address 
basic mechanisms of the charged particle acceleration, the collisionless shock wave and magnetic 
reconnection and vortex dynamics properties relevant to the problem of ultrarelativistic particle 
acceleration. 
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I. INTRODUCTION 



High-power laser facilities have made unprecedented progress in recent years and the near- 
est future their radiation may reach intensities of 10 24 W/cm 2 and higher [1]. As a result of 
laser technology progress the laser-matter interaction entered regimes of interest for astro- 
physics. Typically in the course of laser irradiation of targets shock waves are generated; the 
target compression is accompanied by the Rayleigh- Taylor (RT) and Richtmayer-Meshkov 
(RM) instability development; collimated plasma jets are observed; the matter equation 
of state (EOS) acquires new properties under extreme pressure, density and temperature 
conditions; the laser plasma emits high energy charged particle beams and high- and low- 
frequency powerful electromagnetic radiation. Gathering of these facts principal for both 
space and laboratory physics has initiated works in the so-called laboratory astrophysics [2] 
with the aim to model the processes of key importance for the space objects under laboratory 
conditions. Concerning the laser facilities, the present day laser systems can be subdivided 
into two categories. The first category includes lasers with a relatively long pulse of pico- 
and nanosecond duration and generally low repetition rate. These high energy and power 
laser facilities have been mainly developed for purposes of inertial confinement fusion with 
the laser pulse and target parameters corresponding to the collisional hydrodynamics phe- 
nomena [3]. In context of laboratory astrophysics they are used for experiments on shock 
waves, including the radiative shocks and RT&RM instability, the jet formation, and the 
EOS studies. The second category includes table top size lasers, whose pulse duration is of 
the order of a few tens of a femto-second with high repetition rate [4]. Due to ultra short 
pulse duration and high contrast, these relatively moderate energy lasers can produce ex- 
tremely high power and relativistically high intensity electromagnetic pulses. However, the 
role of both kinds of laser systems is complementary for the development of experimental 
facilities for the purposes of relativistic laboratory astrophysics. 

Generic questions for astrophysics such as whether we are living in the Universe or in 
the Multiverse [5], related discussions of the inflation era in the Multiverse evolution [6] and 
probing our world's dimensions are related to quantum gravitation physics and deal with the 
observational cosmology, in particular with an analysis of the cosmic black body radiation, 
the nuclear synthesis and Type I supernovae radiation (see [7]), are yet out of the energy 
range accessible with present day lasers. The quantum gravitation energy scale is given by 
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the Planck energy, ^hc 3 /G m 10 19 GeV, which corresponds to the mass ^/hc/G m 10~ 5 g 
and the length ^/ftG/c 3 xs 10~ 33 cm. In quantum field theory the unification energy scale 
corresponds to 10 16 GeV [8, 9]. These energy frontiers are yet well above of nowadays laser 
pulse energies. Fortunately, new physics such as the Higgs boson detection and exploration 
of the physics beyond the Standard Model is anticipated to be met at a substantially lower 
energy level in the range of several TeV in the experiments planned with the Large Hadron 
Collider (LHC), as summarized in Ref. [9]. If relativistic laser plasmas can provide the 
charged particle acceleration up to the TeV energy level, laser accelerators will make a 
considerable impact to high energy physics, to finding answers on black hole and brane 
production under the terrestrial conditions [10], to test causality [11] and to study the 
quark-gluon plasmas [12]. 

We may see that the main field of studies of astrophysical phenomena with high power 
lasers lies in the electrodynamics of continuous media in the relativistic regime [13]. Since 
matter irradiated by ultrastrong electromagnetic waves (EMW) is ionized during a time 
interval comparable with the wave period and becomes a plasma and under astrophysical 
conditions approximately 95% of barionic part of matter is in the plasma state, the object 
of our studies is the relativistic laser and astrophysical plasma. 

If we address to the problems of contemporary relativistic astrophysics, first of all ques- 
tions on the mechanisms of the cosmic ray acceleration and on the properties of strong 
EMW interaction with relativistic plasmas attract our attention [14]. In space plasmas 
basic mechanisms of charged particle acceleration are connected with the reconnection of 
magnetic field lines, which is accompanied by the strong and regular electric field genera- 
tion (it occurs in the planet magnetosheres, in binary stellar systems, in accretion disks, 
in the magnetar magnetospheres, etc.) and with collisionless shock waves, at the fronts of 
which the charged particle acceleration occurs (this happens in interplanetary space, during 
supernova explosions, in colliding galaxies, etc.) [14, 15]. 

The laser accelerator development relies upon the fact that under the terrestrial labo- 
ratory conditions presently one of the most powerful sources of coherent electromagnetic 
radiation is provided by lasers [13]. Wakefield accelerators, [16] and [17], presently provide 
the most advanced schemes for electron acceleration and they may be suggested to be good 
candidates for the charged particle acceleration in space [14, 18]. One of the efficient mech- 
anisms of ion acceleration in laser plasmas utilizes the radiation pressure of electromagnetic 
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waves interacting with plasmas (see Refs. [19] and [20] ). Radiation pressure is a very effec- 
tive mechanism of momentum transfer to charged particles. This mechanism was introduced 
long ago [21] and physical conditions of interest range from stellar structures and radiation 
generated winds (see e.g. Refs. [22]), to the formation of "photon bubbles" in very hot 
stars and accretion disks [23], to particle acceleration in the laboratory [19, 24, 25], see in 
addition Refs. [26, 27], and in high energy astrophysical environments [28]. 

Utilization of the plasma nonlinear properties for the electromagnetic wave intensification 
can result in much higher intensity and power. In this case a fundamental role is played 
by relativistic mirrors, which are thin electron sheets induced by the laser radiation moving 
with a speed close to the speed of light in vacuum, as proposed in Ref. [29]. We note 
the fruitfulness of the relativistic mirror concept for solving a wide range of problems in 
modern theoretical physics. Relativistic mirrors are important elements in the theory of the 
dynamical Casimir effect [30], with regard to the Unruh radiation [31] and other nonlinear 
vacuum phenomena [32-34, 36, 37]. Relativistic mirrors made by wake waves may lead to an 
electromagnetic wave intensification resulting in an increase of pulse power up to the level 
when the electric field of the wave reaches the Schwinger limit [38] when electron-positron 
pairs are created from the vacuum and the vacuum refractive index becomes nonlinearly 
dependent on the electromagnetic field strength. In quantum field theory particle creation 
from the vacuum attracts a great attention, because it provides a typical example of non 
perturbative processes [39] . Nonlinear QED vacuum properties can in future be probed with 
such strong and powerful electromagnetic pulses. 

If we trace a relationship between astrophysics and laser physics, we can see a num- 
ber of publications devoted to the laboratory modeling of astrophysical processes [2]. As 
known there has been an interest in modeling space physics with laboratory experiments 
for many years. The first modeling of processes fundamental for space physics in terrestrial 
laboratories has been done by Kristian Birkeland, who more than 100 years ago conducted 
first experiments on studying the auroral regions in the earth magnetosphere [40]. Lat- 
eron progress has been achieved in the laboratory modeling of various processes [41, 42], 
including the magnetic field reconnection [43], collisionless shock waves [44], which provide 
mechanisms for charged particle acceleration under various astrophysical conditions (see Ref. 
[14])- 

In the present paper we address plasma processes relevant to space physics, which occur 
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in the relativistic and collisionless regimes. 

II. DIMENSIONLESS PARAMETERS THAT CHARACTERIZE THE INTERAC- 
TION REGIMES OF HIGH INTENSITY ELECTROMAGNETIC WAVES WITH 
MATTER 

A. Principle of Qualitative Scaling 

Laboratory experiments for studying astrophysical phenomena are of two types [48] . The 
first type of experiments can be referred to as configuration modeling, which is aiming at 
simulating the actual configuration of a system, e.g. the whole Earth's magnetosphere (for 
example see Ref. [45], where the results of the laser-plasma experiments on the simulation 
of the global impact of the coronal mass ejections onto the Earth's magnetosphere are 
presented). The second type of experiments corresponds to process simulation, i. e. they are 
aiming at studying the properties of physical processes relevant to astrophysical phenomena 
[2] . There are a number of nonlinear plasma physical processes that require their clarification. 

Physical systems obey scaling laws, which can also be presented as similarity rules. In 
the theory of similarity and modeling the key role is played by dimensionless parameters 
that characterize the phenomena under consideration [46]. The principle requirement of 
the laboratory modeling is the equality of the key dimensionless parameters in the modeled 
processes. In cases of modeling astrophysical phenomena where this equality can hardly be 
respected, instead the principle of limited similarity (PLS) or principle of qualitative scaling 
has been formulated in Refs. [47] and [48]. According to the PLS those dimensionless 
parameters, which are relevant in a certain context and which are much larger or smaller 
than unity under astrophysical conditions must retain this property (i.e. be much larger 
or smaller than unity) in the laboratory experiments modeling the astrophysical process. 
Below we present the key dimensionless parameters that characterize the high intensity 
electromagnetic wave (EMW) interaction with matter (see also Refs. [49, 50]). 

B. Parameters of Strong EMW Propagating in Plasmas 
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The intensity of an electromagnetic wave pulse is defined by its electric field amplitude 
through the expression: I = cE^/An, which is related to the Poynting vector 

P=±[ExB]. (1) 

The power of the EMW is equal to the integral over its transverse cross-section S, 

P=-^-j([ExB]-n)dS = IS. (2) 
s 

The time integral of the power gives the pulse energy, S = Vr p , where t p , is the pulse 
duration. Other important parameters are the pulse frequency, uq, which is related to its 
wavelength, Ao = 27rc/c<Jo, and the pulse polarization. 

The first of the dimensionless parameters which characterizes the EMW packet is the 
ratio of the pulse length, l p = cr p , to the radiation wavelength, Ao- We shall denote this 
ratio as N p = I p /\q. It is equal to the number of wavelengths per pulse, and is Lorentz 
invariant. 

If the EMW intensity is relatively low, irradiated matter is not ionized. We notice that 
the typical energy of a photon in the laser parameter range with wavelengths in the micron 
range is of the order of one electron-volt and is substantially smaller than the binding 
energy of an electron inside an atom, huj <C Wf,, i.e. it is smaller than the atomic ionization 
potential. In this case the characteristic dimensionless parameter of the interaction is the 
ratio between the amplitude of the electric field in the laser pulse, E = y/Anl/c, and the 
atomic electric field, E a . The latter is equal to the electric field of the proton at a distance 
of a Bohr radius, as = h 2 /m e e 2 ~ 5.3 x lCT 9 cm, i.e. E a = e/a 2 B = m 2 e b /h i . The electron 
binding energy is W b = h 2 /2a 2 B m e and it corresponds to the frequency u a = W b /h. The 
above condition, Hcuo <C W b , is equivalent to the inequality hu /Wb = w /w a <C 1. The 
dimensionless parameter 

Ep _ h A E 

E~ a ~^ U 
becomes equal to unity for a laser radiation intensity equal to m A e w c/ATrh s ~ 10 16 W/cm 2 . 
For small but finite values of this parameter, i.e. in the limit / < 10 16 W/cm 2 , the atom is 
not ionized, unless the multiphoton processes come into play, the EMW-matter interaction 
can be described within the framework of perturbation theory. When the parameter E /E a 
approaches unity, the potential inside the atom changes its form and the so-called tunnel 
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ionization becomes possible. The tunnel ionization probability is given by the Keldysh 
formula [51] 



w = u a exp -- — f(~f K ) , (4) 
nujQ 

where the function fi^x) ~ 27^/3 for j K <C 1 and fi'Jx) ~ hi 2 r y K — 1/2 for ^ K >> 1. The 
adiabatic parameter ^ K is defined as 



eE y «o m e c 
Here introduced is the EMW dimensionless amplitude, 



7^ = ^o — = \l o o - (5) 



eE 

a = . (6) 

m e ujQC 

In the limit 7^ <C 1, i.e. for a relatively strong electromagnetic wave, Eq. (4) corresponds 
to the ionization probability by a constant electric field, 

W = 2 ^t 6XP ("^)- (7) 
For intensities larger than mfe w c/4irh 8 ~ 10 16 W/cm 2 the deformation of the potential 

inside the atoms caused by the laser pulse field becomes so strong that the electron energy 

level becomes larger than the maximum value of the potential. As a result, the electron 

appears as if in a free state and leaves the atom. Due to the periodicity of the electric field, 

there is a probability that the electron will return after a half of the wave period. Recollisions 

with the ions lead to the generation of high order harmonics [52]. However, for a very strong 

electromagnetic wave the effects of the wave magnetic field decrease this probability. In this 

case the matter becomes ionized in one optical period and plasma processes start to play a 

key role. 

Under the action of the electromagnetic wave the plasma electrons oscillate at the 
wave frequency. In the limit v « c their quiver velocity is approximately equal to 
Ve — eE Q /m e Uo. In the non-relativistic limit, when Ve/c <^ 1 or a 1, the electron 
quiver amplitude is smaller than the laser wavelength, Ao- Under the action of the electro- 
magnetic wave, given by the vector potential A±(x — ct), the electrons oscillate at the wave 
frequency. From the equations of the motion we obtain that the transverse component of 
the generalized momentum p± — eA±(x — ct)/c is constant. The particle energy and the 
longitudinal momentum component are related as [53] 

^jml^+pl+p\-p\\c = h. (8) 




b) , 




FIG. 1: Projectories of the charged particle trajectory, when it interacts with the elliptically 
polarized EMW. 

In the reference frame where the particle was at rest before interaction with the laser pulse, 
the particle kinetic energy 



K 



m e c 



y/1 + Go/m e c) 2 - l) 



(9) 



and momentum p = (p\\,p±) are given by expressions K = m e c 2 \a±\ 2 /2, p± = m e ca±, 
p\\ = m e c\a±\ 2 /2. Here a± = eA±(x — ct)/m e c 2 . For |oj_| > 1 the particle acquires a 
relativistic energy, and the longitudinal component of its momentum is larger than the 
transverse component. Fig. 1 shows a typical trajectory of the charged particle in the 
electromagnetic wave. 

The EMW behavior in a plasma differs from its behavior in vacuum, and depends on 
the electron density. In a plasma with a density n, a displacement of the electrons with 
respect to the ions generates the electric field. Its ratio to the laser electric field is E/Eq = 
Anne 2 /m e uj 2 = (uj pe /oJo) 2 = n/n cr , where uj pe = Anne 2 / m e is the Langmuir frequency and 
n cr = irizUj 2 ,/ Aire 2 is the critical density. The dimensionless parameter 



pe 
UJ 



cr 



1 Aime 2 
m e u}Q 



(10) 



is a measure of the plasma collective response to a periodic electromagnetic field. 

When an EMW propagates through a plasma, its group velocity, v g = du/dk, and phase 
velocity, v p h = u/k, are not equal to each other and are related as v g v p h = c 2 . While in 
vacuum the dispersion equation for the frequency, uj, and wave vector,/;;, takes the form 
uj 2 = k 2 c 2 , in a plasma it becomes u 2 = k 2 c 2 + u 2 e . This dispersion equation can be 
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rewritten as k = ^Juj 2 — ou pe /c, which shows that an EMW with a frequency below the 
Langmuir frequency cannot propagate through the plasma and that the electromagnetic 
field evanescence length in a high density plasma is of the order of the collisionless skin 
depth, d e = c/u pe , i.e. an overdense plasma with the electron density higher than the 
critical density is not transparent. 

The collective response of the plasma, in addition to the transverse electromagnetic mode, 
exhibits longitudinal plasma oscillations, i.e., Langmuir waves. The electric field in a Lang- 
muir wave oscillates with frequency uj = uj pe . The group velocity of the Langmuir waves 
vanishes, v g = du pe /dk = 0, and their phase velocity is determined by the wave number. 

Relativistic effects change the dispersion equation due to the dependence of the Langmuir 
frequency on the wave amplitude. As found in Ref. [54], the frequency of a longitudinal 
wave depends on its amplitude a^ — eE / m e u pe c as uj ~ oj pe (l — 3a|/4) for < 1 and as 
uj ~ uj pe j y/8a,L in the case aj_, >> 1. 

For a circularly polarized electromagnetic wave the dispersion equation takes the form: 

V 1 + a l 

We see that the effective critical density increases as the EMW amplitude grows, i.e., the 
plasma is more transparent to high intensity electromagnetic radiation. 

Large amplitude, finite length pulses of electromagnetic and Langmuir waves do not 
propagate independently since they are coupled by nonlinear processes. The Langmuir wave 
that is generated by an ultra short laser pulse, being left behind in the plasma and thus 
called the wake wave. It is of special interest since the structure of the electric field of this 
wake wave is favorable for charged particle acceleration. In a low density plasma the phase 
velocity of the wake wave can be very close to the speed of light in vacuum. In analogy to 
linear accelerators that use electric fields in the radio-frequency range in Ref. [16] it was 
proposed to use the wake field for charged particle acceleration. 

The dimensionless amplitude, Eq. (6), is equal to the electron quiver momentum nor- 
malized to m e c. For a pulse with an intensity corresponding to ao > 1, relativistic effects 
must be taken into account. The intensity of a linearly polarized electromagnetic wave can 
be written via a as 

/ 1.37X10" xaJxteYi!!. (12) 

2 Aq r e \ *o J cm 2 
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If the wave is focused into a one wavelength spot, this intensity corresponds to the power 
V = <2q x 43 GW. At present laser intensities have reached a level above 10 22 W/cm 2 [55]. 

When the electron energy approaches 3m e c 2 , electron-positron pairs are generated during 
electron-nuclei collisions, [39]. The cross section of this process is given by 



Here Ze is the nucleus electric charge and a = e 2 /hc = 1/137 is the fine-structure 
constant. Positron generation in a plasma has been discussed in a number of publications 
(e.g. see [56]) and was observed in the terawatt laser plasma interaction experiments [58]. 
We note a discussion of the pion and muon production in electron-positron and gamma 
plasmas [57]. 

C. Interaction of EMW with Plasmas in the Radiation-Dominated Regime 

The dimensionless parameters characterizing the electromagnetic emission by an elec- 
tron are the ratio between the classical electron radius and the electromagnetic wavelength, 
r e/Ao = e 2 u;o/2vrm e c 3 , and the ratio between the photon energy and the electron rest mass 
energy, hwo/m e c 2 . 

When an electron moves under the action of the electric and magnetic field of a wave, 
it emits electromagnetic radiation. The intensity of this radiation is given by the formula 
W = (2e 2 /3m 2 c 3 )(cip At /rfr) 2 , where is the particle 4-momentum and r is its proper time. 
When an ultrarelativistic charged particle moves along a circular trajectory in a circularly 
polarized electromagnetic wave, the radiation intensity is W = (47rr e /3Ao) ujq m e c 2 dQ. We 
see that the relative role of the radiation damping force is determined by the dimensionless 
parametere rad , which is equal to 



By comparing the energy radiated by the particle per unit of time with the maximum 
energy gain in the electromagnetic wave dtS = uom e c 2 ao, we obtain that the radiation 

1 /3 

effects become dominant at a > a ra d = e rad , i. e. in the limit / > 10 23 W/cm 2 for 1 /iin 
wavelength laser [20, 59]. In the limit of a relatively low amplitude laser pulse, ao <C a ra d, the 
momentum of an electron moving in a circularly polarized electromagnetic wave in a plasma 




(13) 
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scales with the laser pulse amplitude as p = m e ca , while in the limit ao ^> a ra d, it scales as 
p = m e c(a /e rad ) 1/4 . 

Quantum effects become important, when the energy of the photon generated by Comp- 
ton scattering is of the order of the electron energy, i.e. hu m ~ £ e . An electron with energy 
S e = / -fm e c 2 rotates with frequency ujq in a circularly polarized wave propagating in a plasma 
and emits photons with frequency uj m = 7 3 u;o- We obtain that quantum effects come into 
play when 7 > 7q = ^/m e c 2 /fko . For an electron interacting with one-micron laser light we 
find 7q ~ 600. From the previous analysis we obtain that the quantum limit on the electron 
gamma factor corresponds to 

2e 2 m e c 

The energy flux reemitted by the electron is equal to e(E ■ v ) = e ra d^o / y 2 P±/'me- The total 
scattering cross section defined as the ratio of the reemitted energy to the Poynting vector 
P = cEq/Att, is given by 

' = a *rrh? (16) 

where the Thomson scattering cross section is <tt = 87rr 2 /3 = 6.65 x 10~ 25 cm 2 . We see that, 
as the wave amplitude increases in the range 1 < a < a r ad, the scattering cross section 
increases according to the law a = cj t {\ + a 2 ,) and reaches its maximum a = (y T a^ rad at 
a pa a ra( i; for a ^> a r ad, h decreases according to the law o = aTa^. ad /aQ. In Fig. 2 we show 
the scattering cross section dependence on the EMW amplitude and wavelength. 

In the radiation-dominated regime of the EMW interaction with charged particles, i.e. at 
do > o,rad, the emitted gamma quanta can produce secondary electron-positron pairs, which 
in turn emit gamma ray photons, producing an avalanche of 7 rays and electron-positron 
pairs 



D. Probing Nonlinear Vacuum 

When the amplitude of the electromagnetic wave approaches the critical electric field of 
quantum electrodynamics (also called the "Schwinger field"), vacuum becomes polarized and 
electron-positron pairs are created in vacuum [39, 61]. On a distance equal to the Compton 
length, Ac = h/m e c, the work of the critical field on an electron is equal to the electron rest 
mass energy, m e c 2 , i.e. cEqed^c = m e c 2 . The dimensionless parameter 
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-A- = i*£ (17) 



e 



becomes equal to unity for an electromagnetic wave intensity of the order of 

, '!^,„«10»I (18) 
r e X c 47r cm 2 

For such ultrahigh intensities the effects of nonlinear quantum electrodynamics play a 
key role: an electromagnetic wave excites virtual electron-positron pairs. An observable 
manifestation of this process could be detection of light birefringence during the propagation 
of an electromagnetic wave in a strong electric or magnetic field in vacuum. The cross section 
for the photon-photon interaction in the limit hu <C m e c 2 is given by 

^77-77 ioi25vr 2 6 Vm e cV ' 1 ' 

where ftio is the photon energy (see [39]). This cross section reaches its maximum, 
Cmax ~ lCT 20 cm 2 , for hu ~ m e c 2 , i.e. for the interactions of photons in the gamma range. 
Also attention is focused on the process of electron-positron pair creation in vacuum by 
an electromagnetic wave. For an electric field small compared to Eqed, this process is 
sub-barrier, similarly to the tunnel ionization of atom discussed above [see Eq. (4)]. The 
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FIG. 3: Various regimes of relativistically strong EMW interaction with plasmas. 

probability of electron-positron pair creation per unit volume and per unit time is exponen- 
tially small and is given by 



ac \ ( E \ ( E QED 



Here Ac* = h/m e c is the Compton length and a = e 2 /hc = 1/137 is the fine structure 
constant. 

We may formally estimate the number of electron-positron pairs produced by a 10 fs long 
laser pulse in a volume V = A 3 = 10 _12 cm 3 as N± = wVt p . It is easy to show that N± is 
equal to one pair for a laser intensity equal to I = 10 26 W/cm 2 (a more detailed description 
of this process can be found in [65] and in Refs. [39, 62, 63]). Obviously, this latter number 
is overestimated because the minimum needed energy is by many orders of magnitude larger 
than the total energy of the laser pulse. At intensities of the order I = 10 30 W/cm 2 Eq. 
(20) is not applicable and a depletion of the laser pulse must be taken into account. The 
electromagnetic pulse depletion due to its energy conversion into electron-positron pairs has 
been studied in Ref. [64]. 

The nonlinear dependence of the vacuum susceptibilities on the electromagnetic-field 
amplitude results in the finite value of the Kerr constant of vacuum. It can be found to be 

Kk = — A" (21) 
907r m e cA 

The Kerr constant in vacuum for Ao = lfim is of the order of 10 2 cm 2 /erg, which is a factor 
10 20 smaller than for water. As shown in Ref. [32], in a QED nonlinear vacuum two counter- 
propagating electromagnetic waves mutually focus each other. A nonlinear modification of 
the refraction index in vacuum within the framework of the Heisenberg-Euler approximation 
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is characterized by the critical value of the electromagnetic wave power 

= (22) 

When the electromagnetic wave power exceeds this value, the cross modulation nonlinear 
effects affect the wave propagation. We see that the critical power, V cr , depends only on 
the laser pulse wavelength, A , and on fundamental constants. It is easy to show that for 
A = 1/irn the critical power V cr = cE 2 w 2 /A, where w is the laser beam waist. For the 
mutual self-focusing V cr = 2.5 x 10 24 W can be found to be for A = 1/im. 

Nonlinear modifications of the vacuum refraction index lead to the vacuum birefringence 
[32], to the four- wave interaction [33], to the high order harmonic generation [34], and to the 
laser-photon splitting and merging [35] (see also review articles [13], [36, 37]). According 
to Ref. [31] the Unruh radiation intensity of the electron moving in the field of a strong 
electromagnetic wave becomes comparable with the nonlinear Thomson scattering intensity 
under the condition 4ira ftw /m e c 2 ~ 1. The multi-photon Compton scattering during the 
collision of counter-propagating laser beams and ultrarelativistic electron bunches leading 
to the gamma quanta generation 

e~ + nhjjQ — > hu-f, (23) 

with their subsequent interaction with the laser light accompanied by the electron -positron 
pair creation in vacuum via the Breit- Wheeler process 

hu]~ f + nhuj — > e~ + e + (24) 

has been investigated in Ref. [66]. In Ref. [67] the cross section of the Breit- Wheeler process 
on the laser pulse intensity has been investigated. 

Various regimes of the relativistically strong EMW interaction with plasmas are illus- 
trated in Fig. 3. 

E. EMW Parameters under Space Plasma Conditions 

In one of the first works on the charged particle acceleration by strong electromagnetic 
waves in astrophysical plasmas, pulsars [68] have been considered as sources of ultraintense 
radiation [70]. Pulsars are considered to be oblique rotators with non-parallel rotation and 
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FIG. 4: Pulsar magnetosphere. The inset: The Crab pulsar [71]. 



magnetic dipole illustrated in Fig. 3. The power of magneto-dipole radiation is 

given by the expression 



W 



2/i 2 sin 2 x^p 
3^ ■' 



(25) 



where \x is the magnetic momentum, \ is an angle between the rotational and magnetic 
dipole axes, and Qp is the pulsar rotation frequency. Even for parallel magnetic and angular 
moments, i.e. for % = 0, the expression W = (2/3)/i 2 f2p/c 3 gives the pulsar electromagnetic 
energy losses, as it follows from the theoretical model of the pulsar magnetosphere [69] . The 
magnetic moment is related to the pulsar magnetic field and radius as fi ~ Br p. For typical 
values of rp = 10 6 cm and B = 10 12 G we obtain \x = 10 30 G cm 3 . The electromagnetic wave 
intensity at the distance r is equal to I = W/^nr 2 . At the wave zone boundary, r = c/flp, 
the dimensionless amplitude of the electric field is 



ap 



(26) 



For the Crab pulsar with the rotation frequency Qp = 200 s" 1 we find ap = 2 x 10 10 . 

According to Ref. [20, 59] in the limit of high radiation intensity the effects of the radia- 
tion damping should be incorporated into the theory of the electromagnetic wave interaction 
with plasmas. A dimensionless parameter, 



^■rad 



2e 2 Q P 
3m P c 3 



(27) 



— 1/3 

gives a value of the wave amplitude, a rad = e ra j , above which the radiation damping cannot 
be neglected. For Qp = 200 s _1 this yields a ra d = 10 7 , which is substantially less than above 
found value of ap = 2 x 10 10 . 
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In the case of laser - plasma interaction for a typical laser wavelength of l//m the dimen- 
sionless amplitude a rad corresponds to an intensity of the order of 10 23 W/cm 2 which can 
be achieved by tight focusing of the PW power laser beams onto the one-lambda size focus 
spot. We see that the laser plasmas can be used for modeling the radiation damping effects, 
which are important for relativistic astrophysics. 

III. ACCELERATION OF CHARGED PARTICLES IN THE EMW INTERAC- 
TION WITH PLASMAS 

General requirements for the laser accelerator parameters are principally the same as 
for standard accelerators of charged particles [73], i. e. they should have a reasonable 
acceleration scale length, a high enough efficiency and the required maximal energy, a high 
quality, emittance and luminosity of charged particle beams. In the 1940-s Enrico Fermi 
paid attention to the high energy limit of « 1 PeV= 10 15 eV for accelerated particles, which 
could be reached under terrestrial conditions, when the accelerator size is limited by the 
equator circumference. These limitations resulted in the 1950-ties in the proposal to use 
collective electric fields excited in a plasma (collective methods of acceleration) in order to 
accelerate charged particles [24]. 

A. Electron Accelerator 

Wakefield acceleration has been proposed in Ref. [18] for the generation of ultra high 
energy cosmic rays. Below we describe the wake field acceleration mechanism using as an 
example the LWFA scheme. 

Under the condition of minimum laser energy the one stage LWFA accelerator scaling is 
described as it follows [16]. The electric field in a plasma has the form of a wave propagating 
with a phase velocity, v p h,w- A gamma factor corresponding to the wave phase velocity is 
given by the expression 'jp^w — (1 ~ v ph,w I ,( ?)~ X ^ 1 ''■ A condition of the wake wave syn- 
chronization with the driver laser pulse yields v p h,w = v 9t i as , where v g j as ps c(1 — 00^/200%) 
is the laser pulse group velocity. The wavelength of the weakly nonlinear wake wave is 
\ = Ao7p/ l ,w- Assuming the electrostatic potential in the wake is equal to m e c 2 /e, we 
obtain for the fast electron gamma factor 7 e = 2ji hW . The acceleration length is given 
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FIG. 5: Schematic view of the standard linear accelerator of charged particles, a) and the LWFA, 
b). 

by l acc = \lph,wi i- e - lace = ^o1ph,w This gives a relationship between the acceleration 
length and the fast electron energy: l acc = Ao7e //2 . For Ao = 1 jttm and 7 e = 10 6 we obtain 
l acc « 1km [75]. 

In the opposite limit, when the laser transverse width ri as < X p , we need to take into 
account the formation of an electron density cavity moving with the group velocity of the 
laser pulse (see Fig.6, where the wake wave left behind the ultra short laser pulse in the 
underdense plasma is shown). The cavity's transverse size is determined by the laser pulse 
width and its length is of the order of the Langmuir wave wavelength. In this limit, the 
wavelength depends on the amplitude of the Langmuir wave, which in turn depends on the 
laser pulse intensity. For a given laser pulse width the electrostatic potential in the cavity 
is of the order of <fi ~ irnoerf as , and the group velocity of a narrow laser pulse is determined 
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FIG. 6: 2D PIC simulations show that the electrons pushed away by the ponderomotive pressure 
of the laser pulse form the " bow wave" [49] . The electron density distribution (a) clearly shows the 
'swallow-tail' formation during the wake wave breaking in the first period of the wave behind the 
laser pulse. The wakefield (the x-component of the electric field) is excited by the laser pulse in 
an under dense plasma (b). Inset: The bow wave formed by colliding galaxies in the Bullet Cluster 
[72]. 



by its width, i.e. jp^w ~ r ias/\. As a result we find the electron energy scaling [49]: 



7e 



las 
A Q A p 



(28) 



It does not depend on the laser pulse amplitude provided ao > e<j)/m e c 2 . The laser energy 
depletion length in this limit is given by 

2 



In 



a l 



las 



Ao 



(29) 



i. e. it is by a factor a greater than in the ID case. 

Considering the laser electron accelerator for the applications in the high energy physics, 
we find that its parameters should satisfy several conditions in addition to the requirement 
on the maximum particle energy. Parameters of fundamental importance such as the lumi- 
nosity characterize the number of reactions produced by the particles in colliding beams of 
a collider. The luminosity is given by the expression 

N,N 2 



4na v (T z 



(30) 



where Ni and A*2 are the numbers of particles in each of the beams, a y and a z are the 
transverse size of the beam in the y and z directions, and / is the frequency of the beam 
collisions. A product of the luminosity and the reaction cross section gives the reaction 



20 



rate. We see that the luminosity can be increased by increasing the particle number in a 
bunch, Nj, and/or by increasing the repetition rate, /, or by decreasing the transverse size 
of the bunch, <7j, by focusing the particle beam into the minimum size focal spot. The focal 
spot size depends on the beam emittance, which is defined as the surface occupied by the 
bunch in the phase plane ((y,p y ) or (z,p z )). A calculation under the assumption of a round 
transverse shape of the beam (a y — a z — r) results in the expression given by the integral 

e± = — I drdr', (31) 



7T _ 

where r is the transverse size of the bunch and r' = dr/dx = dr/cdt [73]. 

The transverse dynamics of the electron in the field of the wake wave is described by the 
equation (see for example Ref. [76]) 

dr\ 2 



M\^Jt) + < r ^ (32) 
where the electron gamma factor depends on time as 7||(i) = 7 e (l — t 2 /t 2 cc ) with 7 e = ^ 2 hW 
and j p h,w — uo/u pe . In the limit 7|| 1 the electron transverse oscillations are described 
by the dependence of the radial displacement on time: 

( r \ 1/4 

r{t)=rm3 vm) cos 



t 

uj b {t')dt' 



(33) 



where u>b(t) = u;p e 7||(t) 2 is the betatron oscillation frequency and r in j and 7j n j ~ j p h,w 
are the radial coordinate and the electron energy at the injection time, ti n j normalized on 
m e c 2 . Calculating the transverse emittance, we find e± = ni^^Upe/uo) 3 mm mrad, with 
k = fi n j/Xp. The normalized emittance, En = £_i_7e, is equal to e± = nn 2 {u) pe /u)o) mm mrad. 

The electron motion in the electric field of the wake plasma wave is characterized by the 
structure of the phase plane (p x ,X = x — v ph t). A calculation of the energy spectrum of 
fast electrons is done in Refs. [76, 77]). It uses the property of electrons injected at the 
breaking point to move along the separatrix. The electrons, whose trajectories lie on the 
separatrix, where they are uniformly distributed, near the top of the separatrix have an 
electron momentum dependence on the coordinate 

p x = p m (l - X 2 u 2 pe /c 2 a ) = p m (l - t 2 /t 2 acc ). (34) 

The distribution function of the electrons at the target has the form 

f(t,£) = (n b uj pe /V2ca )5(S - £ m (l - t 2 /t 2 acc )). (35) 
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Here S(z) is the Dirac delta function and we have assumed that the electrons are ultrarela- 
tivistic with £ = p x c and £ m = p m c . In order to find the energy spectrum of the electrons 
on the target, we must integrate the function on time in the limits between — t acc and t acc . 
We obtain 



tacc 

dAf^) = rp^_ [ s ( £ _ £m ( 1 _t_X) dt= _ n «±_ 

tacc 



(36) 



2y/2cao y/S m (S m - £) 

i.e. the particle spectrum has a typical form oc l/y/S m — £ near maximum energy. 

Using the above given relationships and estimating a maximum number of particles in a 
bunch as N K 2 n e \p, we obtain the luminosity to be equal to 

£ = 1 ° 34 (io]^) (oX7~) (w) c^V (37) 

Here we assume a round transverse shape of the bunch with r ~ r in j{^ in j/^ e ) 1 ^. Utilization 
of flat bunches with a y ^> a z allows and to achieve larger luminosity [78]. In addition, in 
the case of flat beams the space charge effects and beamsstahlung can be weakened. We 
notice that the radiation damping effects on the LWFA operation have been considered in 
Ref. [79]. 



B. Ion Accelerator 



The mechanism of laser acceleration of ions (protons and other ions) is determined by the 
electric field set up by the space charge separation of hot or energetic electrons and the ions. 
The exact mechanisms entering into the energy transfer from the fast electron to the ion 
energy depends on the specific conditions of the laser-target interaction (see review articles 
[13] and [80]). The proton generation is a direct consequence of the electron acceleration. 

The typical energy spectrum of laser accelerated particles observed both in experiments 
and in computer simulations can be approximated by a quasi-thermal distribution with a 
cut-off at a maximum energy. On the other hand, the applications require high quality 
proton beams, i.e. beams with sufficiently small energy spread A£i/£i <^ lj. For example, 
for hadron therapy it is highly desirable to have a proton beam with A£j/£j < 2% in order 
to provide the conditions for a high irradiation dose being delivered to the tumor, while 
sparing neighboring tissues. In Ref. [81] it was shown that such a required beam of laser 
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accelerated ions can be obtained using a double layer target. Extensive computer simulations 
of this target were performed in Ref. [82] and the results of experimental studies of this ion 
acceleration mechanism are presented in Ref. [83]. 



1. Ion acceleration during plasma expansion into vacuum 

Ion acceleration during the collisionless plasma expansion into vacuum appears to be 
one of the most obvious mechanisms of the ion acceleration [84]. In particular, it has been 
considered as one of the possible acceleration mechanisms in space plasmas [14] . When the 
electrons that have been heated and tend to expand overtake the ions in a relatively small 
volume, the electric neutrality of the plasma breaks and the generated electric field induces 
the ion motion. Although a velocity of the bulk ion and electron motion is of the order 
of the ion acoustic speed, v s = ^T e /m p , a small fraction at the plasma front gains energy 
efficiently Here T e is the electron temperature and m, p is the ion (proton) mass. Under 
the most favorable conditions the ions achieve a kinetic energy which corresponds to an ion 
velocity of the order of the electron thermal velocity, i.e. the maximum ion energy can be 
of the order of m p £ e /m e . We notice here that for the electron distributions with the energy 
cut-off this conclusion requires careful analysis (see Refs. [85, 86]) 

In the limit when the electron energy is relativistic, in order to analyze the ion motion 
one should use the equations of relativistic hydrodynamics, v a T^ = with being the 
energy-momentum tensor, 



Here is the four-dimensional velocity vector, P is the pressure, n is the density in the 
proper frame of reference, W = P + e is the enthalpy with e being the internal energy 
density. 

The self-similar plasma motion depending on the variable x = x /t is described by a 
system of ordinary differential equations 



<9 M (raw") = 0, 



(38) 




(39) 




- 1 



(40) 




\u' - (u X - l)P' = 0, 



(41) 
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with a prime denoting a differentiation with respect to x an d u = v/c. Here we use the 
relativistic equation of state of an ideal gas, 

2 K 3 (m e c 2 /T e ) 

yv = c — — — , (42) 

K 2 {m e c 2 /T e ) 

P = nT e , (43) 
where K n (x) are the modified Bessel functions. In the case of T e =constant we find 

u = 44 

c + c s 

/ \ c/c s 

n = n ^ (45) 



\ C + C s 

with c s being the relativistic speed of sound, 



T e K 3 (m e c 2 /T e ) 
m e K 2 (m e c 2 /T e )' 

In the ultra-relativistic limit the energy spectrum of fast ions has a power-law form, 

dAfM) a e -*»,* r (47) 



2. Radiation pressure dominated regime of the ion acceleration 

A regime of ion acceleration that exhibits very favorable properties has been identi- 
fied in Ref. [19]. Among the wide variety of ion acceleration mechanisms realized in the 
laser-plasma interaction, the radiation pressure dominated ion acceleration (RPDA) has the 
highest efficiency In the RPDA ion accelerator the laser pulse radiation pressure pushes 
forward the irradiated region of a thin foil as a whole. In the relativistic limit, when the 
electrons and ions move together with the same velocity due to a smallness of the electron 
to ion mass ratio, the ion kinetic energy is by a factor rrii/m e times higher than the electron 
energy. In this case the laser pulse interacts with an accelerated foil like with a relativis- 
tic co-propagating mirror. The electromagnetic radiation reflected back by the relativistic 
mirror has almost negligible energy compared to the energy in the incident laser pulse, i.e. 
the laser energy is almost completely transformed into the energy of fast ions. In Fig. 7 
we show results of 3D PIC simulations of this ion acceleration regime. In the course of the 
interaction with a thin overdense plasma slab the multi-petawatt laser pulse forms a cocoon 
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FIG. 7: Results of 3D PIC simulations of the PPDA ion acceleration regime, a) The electromagnetic 
pulse forms a cocoon confining the EMW energy. The right inset shows a cocoon seen in the plasma 
density and an EMW distribution obtained with the 2D PIC simulation. In the left inset we see a 
cocoon formed by the Black Widow pulsar (Pic. NASA), b) Quasi-monoenergetic ion spectrum. 

confining the EMW energy, thus increasing the coupling of the electromagnetic wave with 
the target (see frame a) and the 2D inset). The ions accelerated beyond the GeV energy 
level have a quasi-monoenergetic spectrum (Fig. 7b). We notice that a combination of the 
RPDA mechanism with the use of double layer targets can substantially increase the ion 
acceleration efficiency as demonstrated in Ref. [26]. 

The equations of the irradiated foil motion can be cast into the form [87]: 

| = ^„ (48) 

where pi is a momentum of the foil element, doi is a vector normal to the foil, the index 
i = 1,2,3, and V is the relativistically invariant pressure. In the frame of reference co- 
moving with the foil the radiation pressure is equal to V = E m /2tt, with Em being the 
EMW amplitude. In the laboratory frame of reference we have Eq = E M (ujq/ujm) 2 , where 
ujq and ujm are the wave frequency in the laboratory and boosted frames. They are related 
to each other as ujm/^g = a/(1 — /3) / (1 + /3). Introducing the Lagrange variables r] and 
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£, related to the Euler coordinates as x = x(r],^,t), y = y(r),£,t), z = z(r),£,t), we find 
that the vector normal to the foil surface element is given by dai = e^dx^dx^. Here dxj 
are the vectors directed along the i-axes, e^k is the fully antisymmetric unity tensor, and 
a summation over repeated indices is assumed. Using these relationships we can find the 
equations of foil motion 

dt ~ v £t3k d V d£ ' 1 ' 

9Xi Pi /rrf\\ 

= c i o 3 , =■ ( 50 ) 
at y/m^c 2 + p k p k 

Here u = no/o is the initial surface density, pi = (p x ,P y ,Pz) is the momentum, Xi = (x,y,z) 
is the foil element coordinate and, the index, m p is the ion mass. In the nonrelativistic limit 
for constant pressure V, this system is reduced to the equations obtained in Ref. [88]. 

When a planar foil is irradiated by a normally incident EM pulse, the ions achieve the 
energy 

2w 2 , x 

h l + 2iu' V ; 

where w is the normalized fluence, 



w 



y_ 00 2im lmiC 



In the limit w ^> 1 the resulting ion energy is equal to the ratio of the laser pulse energy, 
£ias i to the total number of accelerated ions, Ntoti ^ ~ £i as /miC 2 N tot . As an example, 
we consider a solid density foil, n = 10 24 cm -3 , of 1 /im thickness irradiated by a laser pulse 
with a transverse size of 100 /im. For the laser pulse energy of the order of 200 kJ we find 
that the accelerated ion energy is equal to 1 TeV with a total ion number of 10 12 . The ion 
acceleration length in this case is approximately equal to l acc ~ 0.5 km. 

In order to achieve high values of the ion bunch luminosity it is highly desirable to 
decrease the transverse bunch size. This can be achieved by modulating the density inside 
the foil, e.g. by a properly modulated laser pulse. The analysis of the linearized equations 
of the foil motion demonstrates the exponential growth of the modulations 

1/3 

— iqi] — ir£ 



x\ 1 \a, (3, t) oc exp 
where 
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(53) 



Trt = tf, — • (54) 

° 6(g 2 + r 2 ) 3/2 A 2 
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and [87] 
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This opens a way for focusing the accelerated ions onto a narrow spot with lower limit given 
by the foil thickness. Using these results we can estimate the RPDA accelerated ion bunch 
luminosity as 

A first indication of the RPDA - regime has been obtained in the experiments [89], 
when a thin foil target has been irradiated by a laser pulse with an intensity approaching 
10 20 W/cm 2 . 



IV. MINI-BLACK-HOLES ON EARTH 



We may see that when LWFA and RPDA accelerators will reach 100 GeV and TeV 
particle energies, which corresponds to the energy range of interest for high energy physics, 
laser accelerators may be considered as a source of ultrarelativistic particle beams with 
parameters comparable to those, which are produced by standard accelerators. 

As an example for the problems in the field of high energy physics and astrophysics which 
may be explored with laser accelerators of charged particles, we note the mini-black-hole 
detection. In the general relativity theory black holes play a fundamental role. The Einstein 
equation [53, 90], 

z rn p 

where m 2 = 1/G is the square of the Planck mass, and is the energy-momentum tensor 
(the units h — c—1 are used), has a Schwarzschild solution for the interval: 

ds 2 = g^dx^dx" = -r](r)dt 2 + dr 2 /r](r) + r 2 dtt 2 (57) 

with 

r](r) = 1 - (2/m 2 p )M/r. (58) 

Here M and f2 are the object mass and the surface element in the 3D space. The metric given 
by this interval has a singularity at r equal to the Schwarzschild radius, R B h = 2M/m 2 . For 
an object with a mass of the order of the solar mass the black hole radius is equal to 2 km. 
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As it was noted above, a black hole with the size of about the Planckian length, 10~ 33 cm, 
has the mass m p = 10~ 5 g, which corresponds to an energy approximately equal to 10 19 GeV. 

The situation may change, if our world's dimension is higher than 3. In accordance with 
modern quantum field theory [91], our world may have higher dimensions (d + 3). The 
additional dimensions are compactified in a sufficiently small scale, R comp - Gravitational 
interaction is present in the whole space due to its universal character. At the small scale 
for r <C R C omp, the gravitational potential of the field produced by an object with mass M 
behaves as 0(r) = Mj +2 M/r^ l+d \ A constant Mf characterizes the gravitational interaction 
in the small scale limit. In the limit of large scale compared with R comp , i.e. for r ^> R comp 
we have the expression 0(r) = M/m 2 r = M/rR d omp Mj +2 . It yields a relationship between 
m P , Rcomp, and M/ being m 2 = M d+2 R comp - The solution of the Einstein equation in space 
with d+3 dimension gives for the interval the formula 

ds 2 = -r](r)dt 2 + dr 2 /r](r) + r 2 dtt 2 d+3 , (59) 

where fld+3 is the surface element and the metric element is equal to 

V (r) = 1 - (R BH /r) d+1 . (60) 

This results in the expression for the black hole radius 

Kj£ = (2/(d + l))MJ (d+1) M/M f . (61) 

For the constant Mf of the order of 1 TeV the black hole radius is equal to Rbh ~ 10~ 4 fm 
(here 1 fm=10~ 13 cm). A probability of a black hole creation is proportional to the cross 
section of this process [10] 

a BH = ttR 2 bh . (62) 

In the collision of 7 TeV proton bunches with a luminosity of the order of the LHC lumi- 
nosity, £ = 10 34 cm _2 s _1 , it is expected that approximately 10 9 mini-black-holes may be 
detected per year. The created mini-black-holes are to be detected by their emission of 
electromagnetic radiation and of elementary particles according to the Hawking mechanism. 
At the end of its evolution the black hole is thought be strings. 

The TeV range laser accelerator of charged particles can generate 10 6 mini-black-hole per 
year, when its repetition rate is 1 Hz. 
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V. FLYING MIRROR CONCEPT OF THE ELECTROMAGNETIC WAVE IN- 
TENSIFICATION 



An electromagnetic wave reflected off a moving mirror undergoes frequency multiplication 
and corresponding increase in the electric field magnitude. The multiplication factor (f + 
/3m)/(1 — /3m) is approximately proportional to the square of the Lorentz factor of the 
mirror, 7^ = 1/ a/1 — [3 2 M) making this effect an attractive basis for a source of powerful 
high-frequency radiation. Several ways have been suggested to extremely high intensity (see 
articles, Refs. [29, 92, 93], [13, 110] and literature quoted in). A specular reflection by a 
sufficiently dense relativistic electron cloud as suggested in Refs. [94] . The reflection at the 
moving ionization fronts was studied in Refs. [95]. 

Here we consider the "flying mirror" concept [29]. It uses a fact that at optimal condi- 
tions, the dense shells formed in the electron density in a strongly nonlinear plasma wake, 
generated by a short laser pulse, reflect a portion of a counter-propagating laser pulse. In 
the wake wave generated by the ultrashort laser pulse electron density modulations take the 
form of a paraballoid moving with the phase velocity close to the speed of light in vacuum 
[98]. At the wave breaking threshold the electron density in the nonlinear wake wave tends 
towards infinity. The formation of peaked electron density maxima breaks the geometric 
optics approximation and provides conditions for the reflection of a substantially high num- 
ber of photons of the counterpropagating laser pulse. As a result of the electromagnetic 
wave reflection from such a "relativistic flying mirror", the reflected pulse is compressed in 
the longitudinal direction, which is a consequence of frequency upshifting. The paraboloidal 
form of the mirrors leads to a reflected wave focusing into the spot with the size determined 
by the shortened wavelength of the reflected radiation (see Fig. 8). This mechanism al- 
lows to generate extremely short, femto-, atto-, zepto-second duration pulses of coherent 
electromagnetic radiation with extremely high intensity, which pave the way for studying 
such nonlinear quantum electrodynamics effects as the electron-positron pair creation and 
nonlinear refraction in vacuum. 

The key parameter in the problem of Flying Relativistic Mirror (FRM) is the wake wave 
gamma factor, ^ p h,w- According to the special theory of relativity [96], the frequency of 
the electromagnetic wave reflected from FRM increases by a factor approximately equal to 
^lph,w A number of back reflected photons is proportional to j~^ w (for details see Ref. 
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FIG. 8: Flying Mirror Concept, a) The reflection of EMW at the relativistic mirror results in a 
frequency upshifting and compression of the wave, b) Paraboloidal modulations of the electron 
density in the plasma wake wave, c) The electric field pattern of the laser pulse driver and of the 
reflected EMW. Inset: The reflected electromagnetic pulse frequency is upshifted, it is focused and 
its intensity increases. 

[97]), which results in the reflected light intensification [29] 

Ir/Io « TpVOW 2 , (63) 

where S is the transverse size of the laser pulse incident on the FRM. The reflected pulse 
power increases as V r = Vojp^w- 

Using the expression for the reflected pulse intensity (63), we obtain that the interaction 
of two laser pulses with energies 10 kJ and 30 J, respectively, counterpropagating in a plasma 
with a density « 10 18 cm~ 3 can result in a light intensification of up to ~ 10 28 W cm -2 . This 
corresponds to the generation of an electric field with a value close to the nonlinear quantum 
electrodynamics (QED) limit, Eqed = mlc 3 /eH, when electron-positron pairs can be created 
in vacuum. This QED electric field is also called the " Schwinger field" . 

Experiments utilizing the electromagnetic pulse intensified with the FRM technique may 
allow studying regimes of super-Schwinger fields, when E > Eq ED . This may be possible 
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because the light reflected by the parabaloidal FRM is focused into a focus spot moving 
with a relativistic velocity and is well collimated within an angle « l/%h,w [29]. The 
wave localization within the narrow angle corresponds to the fact that the wave properties 
are close to the plane wave properties to the extent of the smallness of the parameter 
^/lph,w- In this case the second Poincare invariant of the electromagnetic field, B 2 — E 2 , 
has a value of the order of E 2 /^ 2 hW . Therefore the electric field amplitude in the reflected 
electromagnetic wave can exceed the Schwinger limit by 7 p / 1) vk times. We note that a tightly 
focused electromagnetic wave cannot have an amplitude above Eq ED , due to the electron- 
positron pair creation [65] when E — > Eq ED leads to the depletion of the electromagnetic 
wave [64]. 

As it was shown above, the critical power for mutual focusing of two counterpropagating 
EMW is equal to V cr = 2.5 x 10 24 W, which is beyond the reach of existing and planned 
lasers. Fortunately, if we take into account that the radiation reflected by the FRM has a 
shortened wavelength A r = ^o/^p^w an< ^ that ^ s P ower is increased by a factor %h,w, we 
may find that for jp^w — 30, i.e. for a plasma density ~ 3 x 10 17 cm -3 , nonlinear vacuum 
properties can be seen for laser light the incident on the FRM with a power of about 10 
PW. This makes the FRM concept attractive for the purpose of studying nonlinear quantum 
electrodynamics effects. 

Within the framework of the Flying Mirror concept, it has been demonstrated [29] that 
the wavelength of the laser pulse, which has been reflected and focused at the wake plasma 
wave, becomes shorter by a factor and its power increases by a factor 2 r y ph . From this 
it follows that nonlinear QED vacuum polarization effects are expected to be observable for 
50 PW l-/xm lasers. 

A demonstration of the Flying Mirror concept has been accomplished in the experiments 
of Ref. [99]. Two beams of terawatt laser radiation interacted with an underdense plasma 
slab. The first laser pulse excited the nonlinear wake wave in a plasma with parameters 
required for the wave breaking, which has been seen in the quasi-mono-energetic electron 
generation and in the stimulated Raman scattering. The second counter-crossing laser pulse 
has been partially reflected from the relativistic mirrors formed by the wake plasma wave. 
We detected the electromagnetic pulses with a duration of femtoseconds and wavelengths 
from 7 nm to 15 nm. These results demonstrate the feasibility of constructing sources of 
coherent X-ray radiation with the parameters that are tunable in a broad range. 
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VI. RECONNECTION OF MAGNETIC FIELD LINES & VORTEX PATTERNS 



The term " magnetic field line reconnection" refers to a broad range of problems that are of 
interest for space and laboratory plasmas. The results of theoretical and experimental studies 
of magnetic reconnection have been reviewed in many papers and monographs [14, 102, 103]. 
As it concerns relativistic laser plasmas, earlier a conclusion was made in Ref. [104] about 
the important role of the generation of magnetic fields by fast electron currents and their 
reconnection in the relativistic laser-matter interaction regime. The experiments conducted 
in Refs. [105] revealed magnetic reconection phenomena in laser plasmas, when two high 
power laser beams irradiated a thin foil target. 

Processes of reconnection are accompanied by an ultra fast magnetic energy release, 
which is transformed into different forms, such as internal plasma energy, radiation and fast 
particles. 

A. Dimensionless parameters describing the relative roles of nonlinear, dissipative 
and Hall effects 

Inside current sheets, which are basic entities in the reconnection process, as well as in 
the vicinity of shock wave fronts, the effects of dissipation and of nonlinearity play a crucial 
role being comparable in magnitude. Together with the Hall effect, which leads to the 
appearance of small scale structures, these effects violate the freezing of magnetic field in 
plasma motion. 

Magnetic field is frozen in a plasma in the limit of a large Lundquist number S* ^ oo. It 
obeys the equation: 

9,B = vx B, (64) 

which corresponds to the conservation of magnetic field flux through a contour moving with 
the plasma. The dimensionless parameter S is equal to the ratio of two characteristic time 
scales: the magnetic diffusion time r a = l 2 ' jv m and a typical time T\ — 1/va that it takes for 
an Alfven wave to propagate along the distance I; S — r a /rA- Here the magnetic diffusivity 

IS Yfi ^ 2 /47Tcr, where a is the electric conductivity of the plasma, and va = |B|/^/47rp is 

the Alfven wave velocity. 

In the vicinity of the zero point the scale of the field nonuniformity / equals to the distance 
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r from the zero point. The magnetic field and hence the Alfven velocity are proportional to 
r: |B| = hr, va = hr/ y/A.np = fi^r. Here h is a typical value of the gradient of the magnetic 



The measure of the significance of nonlinear effects is given by the ratio SB /hi between the 
magnitude of the magnetic field perturbation SB and the background magnetic field B = hi. 
This ratio depends on the distance from the null point, due to both the nonuniformity of the 
background magnetic field B and to the change of the MHD wave amplitude in the course 
of its propagation. 

If the plasma is pinched by a the quasi-cylindrical electric current I with a radius of the 
order of r, the value of the magnetic field at its boundary is approximately SB = 11 jcr. The 
dimensionless ratio SB /hi is equal to one for r pa r m = yjl/hc. If the electric current has the 
form of a quasi-one-dimensional slab pinch, and if the pinching occurs in the direction of its 
small size, the characteristic value of the magnetic field perturbation is constant: SB = B\\ 
and the ratio e becomes of order unity at the distance = B\\/h. In the approximation 
of small amplitude perturbations, these two types of pinching correspond to the effects of 
the propagation of magnetoacoustic and of Alfven waves, respectively. The magnetoacoustic 
waves focus towards the null line, while the energy of the Alfven waves accumulate near the 
magnetic field separatrices. The values r m and ta determine the size of the region, where 
the magnetoacoustic wave and, respectively, the Alfven one become nonlinear. 

The dimensionless parameters 



determine the relative role of the dissipation and of the nonlinearity effect in the course of 
the current sheet formation due to finite amplitude perturbations of the magnetoacoustic 
and of the Alfven wave type, respectively. 

Now we discuss the relationship between the dimensionless parameter L m and the current 
sheet parameters obtained in the framework of the Sweet - Parker model [106, 107]. In this 
model it is supposed the current sheet has a width b and a thickness a with b ^> a. The 
plasma flows into the current sheet with a velocity v- in ~ v m /a and exits through its narrow 
edges with a velocity v out , which is of the order of the Alfven wave velocity, va ~ &Ab. From 
mass conservation we obtain v m b = av out . From this it follows that the thickness of the 



field. 




(65) 



(r A /r a ) 2 = Bp A /h 2 u m = L A 



(66) 
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current sheet is equal to 

a = y/v m /Q A , (67) 

i.e. of order r a . Estimating the current sheet width as r m , we find that the ratio of its width 
to its thickness is 

b/a = ^JlVt A /hcv m = \flZi- (68) 

Thus, the condition for the formation of a wide current sheet with b ^> a is equivalent to the 
requirement L m ^> 1. Similarly, current sheets are formed in the vicinities of the magnetic 
field separatrices when La 1. 

Considering the case when the Hall effect, i.e. the electron inertia, plays a dominant role 
in the reconnection process, we define the dimensionless parameter which measures the role 
of the Hall effect as a = ah/il A l = c/uo pi l. When the length r H) at which the Hall effect 
starts to be important, is larger than the current sheet thickness, r H /a = ah/ '\Am^A > 
1, the effects of dispersion lead to the formation of small scale structures. In the limit 
r ujb = ah^/chjl /Q,a ^ 1 the pattern of the plasma flow is completely determined by the 
Hall effect. Similar to the way used to define the parameters L m and La, we define the 
dimensionless parameter 

L H = (b/r H ) 2 = E 2 c 2 {l a /h A a 2 v m = E 2 uj 2 Jh 2 v m Vt a . (69) 

When L H 3> 1, nonlinear effects are much stronger than the Hall effect. 



B. Current Sheet 



In a simple 2D configuration the current sheet is formed in the magnetic field described 
by a complex function B(x,y) = B x — iB y = h( of a complex variable ( = x + iy. The 
magnetic field vanishes at the coordinate origin. The magnetic field lines lie on the surfaces 
of constant vector potential, A(x,y) = Ke{h( 2 /2}. They are hyperbolas as we can see in 
Fig. 9 a). This is a typical behaviour of the magnetic field lines in the vicinity of null lines 
(they are the so called X-lines) in magnetic configurations. Under finite time perturbations 
the magnetic X-line evolves to the magnetic configurations of the form B = h(( — b) 1 / 2 , 
which desribes the magnetic field created by thin current sheet between two points ±b [102]. 
The magnetic field lines lie on the constant surfaces of 

A(x, y) = ^Re {(VC 2 ^ 2 ~ Log [( + VC 2 ^} } • (70) 
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FIG. 9: Magnetic field pattern in the vicinity of the X-line (a). Current layer formed in the vicinity 
of the X-line (b). 




FIG. 10: Electric current density distribution inside the current sheet [108]. 

They are shown in Fig. 9 b. The width of the current layer b is determined by the total 
electric current / inside, and by the magnetic field gradient, h. It is equal to 

b = y/UJhc. (71) 

In the strongly nonlinear stage of the magnetic field and plasma evolution a quite complex 
pattern in the MHD flow in the nonadiabatic region near the critical point can be formed, 
with shock waves and current sheets. In Fig. 10 we show the results of the dissipative 
magnetohydrodynamics simulations of the current sheet formation near the X-line. 

C. Magnetic Reconnection in Collisionless Plasmas 

When the Hall effect is dominant, i.e. the electron inertia determines the relationship 
between the electric field and the electric current density carried by the electron component, 
the magnetic field evolution is described by the equation (see [116, 117]) 

d t (B - AB) = V x [(V x B) x (B - AB)] , (72) 
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FIG. 11: Nonlinear stage of the development of the tearing mode instability in a current sheet: a) 
the magnetic field and b) the generalized vorticity distribution at t=0. The same functions at t=8 
in c) and d). 

which corresponds to the condition of generalized vorticity, f2 = B — AB, be frozen into 
the electron component motion with the velocity v e = cV x B/47rnoe. Here the space scale 
is chosen to be equal to the collisionless electron skin-depth, d e = c/u> pe , and the time unit 
is = m e c/eB. The range of frequencies described by the EMHD equations is given by 

10 Bi <U) < U) Be- 
lli the linear approximation Eq. (72) describes the propagation of whistler waves, for 
which the relationship between the wave frequency and the wave vector, is u = |k|(k • 
B )/(l + k 2 ). From this relationship it follows that in a weakly inhomogeneous magnetic 
field the critical points are the points and lines where |B | = or/and (k • B ) = 0. 

The electron inertia effects make the reversed magnetic field configuration unstable 
against tearing modes [118? ], which result in magnetic field line reconnection. The slab 
equilibrium configuration with a magnetic field given by B = B 0z e z + B 0x (y/L)e y , where 
B 0x (y/L) is the function that gives the current sheet magnetic field, is unstable with respect 
to perturbations of the form f(y)exp( , yt + ikx) with kL < 1. For this configuration one 
has (k ■ B ) = at the surface y = 0. The growth rate of the tearing mode instability is 
[117, 119, 120] 7 « (1 - kL) 2 A' 2 /kL 2 . 

In Fig. 11 the results of a numerical solution of Eq. (72) in a 2D geometry with magnetic 
field B(ar, y,t) = (V x a) x + be» are shown. The unperturbed configuration is chosen 
to be a current sheet, infinite in the x-direction, that separates two regions with opposite 
magnetic field. Both the line pattern of generalized vorticity, Q = a — Aa, and of the 
magnetic field show the formation of quasi-one-dimensional singular distributions in the 
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FIG. 12: a) The current sheet break up into two parts separated by the distance 2a(t). b) The 
projections of the trajectory of charged particle accelerating in the vicinity of the masgnetic X-line. 
Inset: The solar flare (YOHKOH image). 

electric current density and in the distribution of the generalized vorticity. The magnetic 
field topology changes, as is seen from Fig. 11. 



D. Charged Particle Acceleration 

A fully developed tearing mode results in a current sheet break up into parts separated by 
a distance 2a, as it is illustrated in Fig. 12a (see. Ref. [102] and literature quoted therein). 
Under the magnetic field line tension the plasma is thrown out. The model magnetic field 



describing this configuration is given by the complex variable function -B(C) = B Q q/ \/a 2 — q 2 . 
The magnetic field lines lie on the surfaces of constant vector potential, 

A(x,y,t) = Re [b ^a 2 (t) - ? 2 } . (73) 

Due to a dependence of the function a on time the electric field parallel to the z axis arrises. 
It is given by 

1 „ 1 Bna(t)a(t) , , 

E(x, y, t) = — d t A = , ° y y . 74 



In the vicinity of the null line we have a quadrupole structure of the magnetic B(() m B^q/a 
field and a locally homogeneous electric field, E w clBq/c. 

The magnetic field reconnection, the study of which has been started by Dungey [121], on 
its initial stage had had as a main goal to explain the generation of suprathermal particles 
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during solar flares and substorms in the earth's magnetosphere. Despite the simplycity 
of the formulation of the problem, it is quite far from a complete solution. Even in the 
test particle approximation, which describes the particle motion in the given magnetic and 
electric fields, the solution of this problem meets serious difficulties [14, 101]. The reason 
of that is due to the fact that in the vicinity of critical points of magnetic configurations 
the standard approximations adopted to describe the plasma dynamics are no longer valid. 
In such regions the drift approximation, i.e., the assumption that the adiabatic invariants 
are constant, can no longer be applied. On the other hand, the particle spends only a finite 
time interval in the nonadiabatic region, since there its motion is unstable. After a finite 
time interval it gets out of the nonadiabatic region, and gets into the drift region as it is 
seen in Fig. 12b. Matching the solution described by the particle trajectories in different 
regions, we can describe the particle motion and hence the acceleration near critical points 
of the magnetic configurations. 

Under the conditions of space plasmas, the radiation losses during the charged particle 
acceleration in the magnetic reconnection processes are caused by the backward Compton 
scattering and by synchrotron losses. A characteristic time of the synchrotron losses for the 
electron with energy S is given by the expression 

TB = (75) 

As it was shown in Ref. [109], during solar flares this effect limits the ultrarelativistic 
electron energy to a value of about several tens of GeV. 



E. Electron Vortices in Collisionless Plasmas 



The vortical fluid motion is well known to be widely present under the earth's and space 
conditions. In laser plasmas, when ultra short and high intensity EMW pulse propagates 
in the collisionless plasmas, it accelerates a copious number of relativistic electrons. The 
electric current of fast electrons produces quasistatic magnetic field, whose evolution results 
in the formation of electron vortex structures. They naturally take a form of the vortex rows 
[111], as it is shown in the LHS inset to Fig. 13. A strong magnetic field in the relativistic 
laser plasma has been detected experimentaly [112]. 

The interacting vortices can be described within the framework of a two-dimensional 
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theoretical model. By taking B to be along the z-axis (B =Be z ), and assuming all the 
quantities to depend on the x, y— coordinates, we obtain from vector equations (72) one 
equation 

d t (AB -B) + {B, (AB -B)} = (76) 
for a scalar function B(x,y,t). Hear 

{f,g} = d x f d y g-d x gd y f (77) 

are the Poisson brackets. Equation (76) is known as the Charney equation [113] or the 
Hasegawa-Mima (HM) [114] equation in the limit of zero drift velocity. In this case linear 
perturbations with the dispersion equation uj = |k|(k-Bo)/(l + /c 2 ) correspond to the Rossby 
waves, the drift waves or to the whistler waves, respectively. 

Equation (76) has a discrete vortex solution, for which the generalized vorticity is localized 
at the points x = x Q : 

n = AB-B = ^/t a <5(x-x Q (t)). (78) 

a 

Solving this equation we find that the magnetic field is a superposition of the magnetic fields 
created at isolated vortices localized at the coordinates x a (t): B = J2 a ^ a 

iT(x,x«(t)) = -^ (|x-x«(t)|). (79) 

Here and below K n (£) are modified Bessel functions. 

The curves x Q (t) are determined by the characteristics of equation (76). The character- 
istic equations have the Hamiltonian form 

•«- T dH 1 J ST {X? Z X jl 

K a X i — Jij — Jij yKaKp , [HO) 

where is the antisymmetric unit matrix. The Hamiltonian is given by 

H = - K a KpK (l a p) /2ir. (81) 

a<!3 

In the case of the Euler hydrodynamics, a point vortex is described by (k q /27t) In |x — 
x Q (t)|, instead of the expression (79) which involves the Bessel function K (|x — x a (t)|). 
The later results in the shiealding of the interaction between vortices at large distances. A 
typical scale length of the problem under consideration in the case of the EMHD vortex 
systems, is equal to the collisionless electron skin-depth, d e = c/u pe . 
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FIG. 13: Instability growth rate on s and q for the vortex row described within the framework 
of the Euler (a) and Hasegawa-Mima (b) approximations. Left inset: The vortex row seen in the 
magnetic field patch distribution in a plasma behind the laser pulse. Right inset: The Mouse pulsar 
(NASA/CXC/SAO, Chandra image of G359.23-0.82 pulsar) [124]. 

Considering the problem of the stability of an infinite vortex chain we assume that all 
vortices have the same absolute intensity and take. In the initial equilibrium the vortices 
have coordinates (Fig. (13)) 

In the case of an antisymmetrical vortex row with a = 1/2, we expect a more complicated 
behavior of the perturbations, compared to that of the symmetrical configuration. As noted 
in Lamb's monograph [115], in standard hydrodynamics the antisymmetrical von Karman's 
vortex row is stable for q/s ~ 0.281, where s and q give a distance between the vortices in 
the unperturbed vortex row along the x and y coordinates. A dependence of the instability 
growth rate on s and q for the vortex row described within the framework of the Euler 
hydrodynamics approximation is shown in Fig. 13a. 

By direct inspection of the row instability described by the Hasegawa-Mima equations 
we can see that for large distance between neighbouring vortices the antisymmetric vortex 
row is stable for 

3s 2 /4 > q > s/2. (82) 

(see Fig. 13b). 
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F. A Role of the Weibel Instability in the Quasistatic and Turbulent Magnetic 
Field Generation 

The quasistatic magnetic field generation in relativistic laser plasmas occurs due to the 
fast electron beam interaction with the background plasma. It can be understood in terms 
of the Weibel instability [122] or in the generic case, in terms of the electromagnetic filamen- 
tation instability. When the fast electron beam propagates in the plasma, its electric current 
is compensated by the current carried by the plasma electrons. A repulsion of the oppositely 
directed electric currents results in the electron beam filamentation and in the generation 
of a strong magnetic field [125]. An electromagnetic filamentation instability leads to the 
generation of a quasistatic magnetic field and is associated with many small-scale current 
filaments [126]. Each filament consists of a direct and of a return electric current which 
repel each other. This produces a strong electric field, which accelerates the ions in the 
radial direction. In the long term evolution, the successive coalescence of the small-scale 
current filaments forms a large scale magnetic structure. This process is accompanied by the 
reconnection of the magnetic field lines, by the formation of current sheets, and by strong 
ion acceleration inside these sheets [127]. 

The filamentation phenomena are of great interest for the explanation of the quasi-static 
magnetic field origin in laser plasmas irradiated by relativistically strong EMW [123]. Coun- 
terstreaming electric current configurations naturally appear in space at the fronts of col- 
liding electron-positron and electron-ion plasma clouds [128] as in the cases of the Galactic 
Gamma Ray Bursts and in shock waves in supernova remnants. The filamentation instabil- 
ity generates the magnetic field required by the theory of the synchrotron afterglow in GRB 
[129]. The Weibel instability has been invoked as a mechanism of the primordial magnetic 
field generation by colliding electron clouds in cosmological plasmas [130]. 

The filamentation instability developing in the vicinity of shock wave fronts together 
with other types of instabilities [131] plays the role of the source of strong electromagnetic 
turbulence invoked in the theoretical models of the Fermi acceleration of cosmic rays [132, 
133]. A realization of the Fermi acceleration mechanism of Type A at the shock wave front 
is discussed below in Section VIII. 
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FIG. 14: Schematic pulsar magneto-sphere, according to Ref. [134]. A rotating relativistic electron 
lump emits electromagnetic radiation by the antenna mechanism. 

VII. RELATIVISTIC ROTATOR 

In Ref. [134] the antenna mechanism of the pulsar radiation emission has been proposed. 
According to this mechanism in the pulsar magnetoshere, which is a rotating magnetic dipole, 
the magnetic dipole interaction with a plasma at the magnetosphere periphery induces strong 
modulations of the electron density, an electron density lump. The phase velocity of the 
electron lump can be arbitrarily close to the speed of light in vacuum. It is directed along a 
circle as illustrated in Fig. ??. As a result of the curvilinear acceleration, the electron lump 
emits a radiation, whose properties are similar to the synchrotron radiation [135]. 

In the context of Relativistic Laboratory Astrophysics it is remarkable that the relativistic 
rotating dipole can naturally be formed in the laser plasma. Laser-plasma interactions 
provide an opportunity to reproduce nonlinear electrodynamics effects under astrophysical 
conditions in the laboratory. In Ref. [136] it is demonstrated that high-power coherent 
synchrotron-like radiation can be generated by the relativistic charge density wave rotating 
self-consistently inside an electromagnetic-dipole solitary wave, dwelling in a laser plasma. 
The relativistically strong laser pulse can generate relativistic EM subcycle solitary waves 
in a plasma [137], as it was indicated by particle-in-cell (PIC) simulations. An analytical 
description of solitons of this type was developed in Refs. [138]. Figure 15 presents the 
structure of electric and magnetic fields inside the soliton [139]. The soliton ressembles 
an oscillating or rotating electric dipole. The toroidal magnetic field, shown in Fig. 15, 
indicates that, besides the strong electrostatic field, the soliton also has an electromagnetic 
field. The electrostatic and electromagnetic components in the soliton are of the same order 
of magnitude. 



42 




- arrui 



FIG. 15: Structure of electric (a) and magnetic (b) fields inside the EM relativistic soliton. Inset: 
The magnetic- and electric-field topology in the TE (with poloidal magnetic field and toroidal 
electric field) and in the TM (with poloidal electric field and toroidal magnetic field) solitons. 

The 3D solitons emit high-frequency EM radiation, whose frequency is much higher than 
the Langmuir frequency [136]. This radiation is emanated from the electron density hump 
rotating in the wall of the soliton cavity, similar to coherent synchrotron-like emission. This 
radiation has the characteristics of a well pronounced outgoing spiral EM wave, Fig. 16 a). 
The emission of the spiral wave correlates to the rotation of the electron density hump in 
the cavity wall, and it leads to the spiral modulations of the electron density (see Fig. 16 
b)). The density hump gyrates in a circle, and the period of revolution is exactly equal 
to the soliton period. The polarization of the spiral wave corresponds to the well known 
synchrotron radiation [135] and the density hump emission is coherent. 

The results of 3D PIC simulations, presented in Ref. [136], distinctly demonstrated 
relativistic rotating dipoles excited by the circularly polarized laser pulse in an underdense 
plasma. The dipoles are associated with the relativistic electromagnetic solitons. 

VIII. SHOCK WAVES 

Phenomena taking place at shock-wave fronts play a key role in various astrophysical 
conditions. The characteristic dimensionless parameters that determine the shock wave 
propagagation are the magnetic Mach number, Ma = vsw/ v a, equal to the ratio of the 
shock wave front velocity, vsw, to the Alfven velocity, v a, the ratio of the gas pressure to 
the magnetic pressure, 8irnT/B 2 , and 9, the angle between the normal to the front and the 
magnetic field. 
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FIG. 16: a) Cross sections of the magnetic field component eBz = m e c in the plane x, y. b) The 
electron density distribution. Right inset: The EM field of the rotating electric charge rotating. 
Left inset: A frequency spectrum of the emitted EM wave. 



FIG. 17: Schematic view of the shock waves and contact dicontinuity in a supernova remnant. In 
the inset: X-ray and optical image of supernova remnant SNR 1987A [148, 149]. 

A. Shock Waves in Supernova Remnants 

The origin of cosmic rays (CR) is one of the most interesting problems in astroparticle 
physics [14, 142]. The observation of ultrahigh-energy cosmic rays indicates that cosmic 
rays exist beyond 10 20 eV and certainly beyond 10 19 eV energies greater than the GZK 
cutoff [143] for the extragalactic sources due to the pionization loss of protons that decay 
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by collision with cosmic microwave background photons. The galactic CR spectra in the 
energy range above a few GeV and below 10 7 GeV are power-laws with the total cosmic 
ray spectrum being 

Jc-llxg- TTpv (83) 
cnrs st GeV 

in the energy range from a few GeV to 100 TeV with k 2.7. Around 10 15 eV (the "knee"), 
the slope steepens from k « 2.7 to k « 3. The energies 10 18 eV correspond to the ultra 
high energy cosmic rays (UHECR), which sources are associated with active galactic nuclei 
(AGNs) [144]. 

For the most advanced theoretical models of galactic cosmic ray acceleration with the en- 
ergy below 10 17 eV the shock waves formed in the supernova explosions are most important. 
This process is related to the nature of collisionless shock waves [131]. 

During explosions of type II supernovae an energy £ tot of the order of 10 51 erg is released. 
The frequency of supernova explosions is about 1/30 per year. Estimates [14] show that 
approximately 2% of the energy of a supernova should be transferred into the cosmic ray 
energy. 

In the initial stage of the evolution of a supernova envelope a system of shocks is formed 
(17). The matter ejected from a star is decelerated and compressed in the inner shock wave. 
Through the circumstellar gas a second shock wave propagates. The matter ejected from 
a star is separated from the circumstellar gas by a contact discontinuity, which is unstable 
with respect to a Rayleigh- Taylor instability. The RT instability leads to the relatively long 
scale modulations of the gas density inside the supernova shells. 

When the mass of the swept interstellar gas becomes larger than the mass ejected from 
the star, the propagation of outer shock in Fig. (17) is described by the Sedov- Taylor self- 
similar solution. The radius of the shock, Rsw, as a function of time is related to the energy, 
Ssn, released in the explosion and to the gas density p by the relation 

Rsw(t) = 1.51 5 t^= 5 -v S wt (84) 

The shock wave velocity, vswif) ~ t~ 3 ^ 5 , decreases with time. At a later time when the 
radiation losses become important, the law of the supernova envelope expansion changes. 
The asymptotic time dependence of the SN envelope radius is given by Rgw(t) ~ t 2 ^ (see 
Ref. [14] and references therein). 
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FIG. 18: Structure of collisionless shock wave front for (3 > (a) and for < (b). 
B. Collisionless Shock Waves 

If the shock wave has a relatively small amplitude, Ma < Mi m 1.5 (the precise value 
depends on (3 and 8), then the front profile is laminar in structure and it is determined by 
a joint action of the dispersion and dissipation on the nonlinear waves propagation. These 
effects are described in the framework of the Korteweg-de Veries- Burgers equation: 



The stationary wave propagating with constant velocity is described by a solution, which 
shows the change of the amplitude of the wave from zero far ahead of the shock wave front, 
to U\ = 2v sw far behind the shock wave front. 

The decay of the oscillation amplitude, with the coefficient equal v/[3, results in the 
decrease of the amplitude of solitons as it is shown in Fig. 18. If dissipation effects are more 
important than the effects of dispersion, vj(3 3> 1, there are no oscillations at the shock 
wave front. More precisely, the decay should be large enough, v 3> v cr , with 



In this case the wave has a monotonous structure. 

In the case of u//3 <C 1, the dispersion effects are dominant and there are many well seen 
solitons near the front. For (3 > the oscillations are localized behind the front (Fig. ??a), 
while for (3 < they are ahead of the front (Fig.??b). 

For example, in the case of the magnetoacoustic shock waves in a plasma, propagating 
almost perpendicularly to the direction of the magnetic field, the dispersion coefficient, 
(3 in v a c 2 /2u)p e , is positive. This means that the oscillations are localized behind the front of 



d t u + ud x u — vd xx u + f3d xxx u = 0. 



(85) 



v, 



cr 



(86) 
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the magnetoacoustic shock wave propagating perpendicularly to the magnetic field. When 
the direction of the magnetoacoustic wave propagation is almost parallel to the direction 
of the magnetic field, the coefficient (3 ~ —v a c 2 /2uj 2 i is negative with uj pi = ^4irne 2 Jrn~i. 
The oscillations at the front of the magnetoacoustic shock wave, propagating quasi-parallely 
with respect to the magnetic field, are localized ahead of the front. 

Dissipation, which determines the distance of the oscillation decay, can be due to anoma- 
lous resistance and viscosity arising from an excitation of the plasma instability, i.e. the 
Weibel instability of counterpenetrating plasmas. If the amplitude of the shock wave is 
large, M A > 3, a high level of turbulent fluctuations of electric and magnetic fields are 
excited ahead and behind the wave front. 

In the laser-plasma physics context, the observation of collisionless shocks was reported 
by several authors [44] , aiming to reproduce astrophysical phenomena in small scale labora- 
tories. However, in general when the shocks were observed with optical probing techniques, 
the front structure could hardly be resolved. In Ref. [145] the propagation in a rarefied 
plasma (n e < 10 15 cm~ 3 ) of collisionless shock waves being excited following the interaction 
of a long (L = 470ps) and intense (/ = 10 15 Wcm~ 2 ) laser pulse with solid targets, has been 
investigated via proton probing techniques [146]. The shocks' structures and related elec- 
tric field distributions were reconstructed with high spatial and temporal resolution. The 
experimental results are described within the framework of the nonlinear wave description 
based on the Korteweg-de Vries-Burgers equation (85). 

C. Diffusive Acceleration of Charged Particles at the Shock Wave Front 

The charged particle interaction with fluctuations of the electric and magnetic field in a 
turbulent plasma may result in particle scattering and diffusion. When the shock wave prop- 
agates in a turbulent medium, an average velocity of electromagnetic fluctuations is different 
in the regions ahead and behind the shock front. Efficiently the particle appears to move 
between semitransparent (due to diffusion) walls with a decreasing distance between them. 
A model transport equation describing the particle convection, diffusion and acceleration 
has a form [14] 

9 t / + div(u/- J DV/) = 
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(Jdiv u - K(pj) f] , (87) 

where f(p,x,t) is the fast particle distribution function, p, x and t the particle momentum, 
coordinate and time, vsw being the speed of the shock wave propagation, and D is the diffu- 
sion coefficient. A term in the right hand side describes regular acceleration or deceleration 
of the charged particles: 

^ = -K(p)-^pdivu. (88) 

The function K(p) corresponds to the Compton and synchrotron losses important for the 
cosmic ray electron component: 

K(p) = -(3 B cp 2 (89) 

with 

Pb = 8 x 1(T 25 + w ph ) (90) 



^8tt p "J eV s 

An average change of the particle momentum proportional to pdivu occurs due to the 
particle bouncing between converging, div u < , or diverging, div u < 0, scattering centres. 

The average particle bouncing between two reflecting plates with distance L as a function 
of time provides a simple example of a dynamic system with conservation of the longitudinal 
adiabatic invariant, J\\ = pL [150, 151]. The phase plane shown in the inset to Fig. 19, 
illustrates the Fermi acceleration mechanism of the first type (type A according to Ref. 
[152]). By virtue of the longitudinal adiabatic invariant conservation, for decreasing distance 
between the plates, dL/dt < 0, the particle momentum grows, i.e. the particle acquires 
energy. 

The velocity distribution in the vicinity of the front of an infinitely thin shock wave, 
propagating from left to right, has the form: u(X) = U\ in the region X > 0, and u(X) = 112 
for X < 0. Here X = x — vswt- The velocities ahead the shock front and behind it are 
related to each other as 

k + 1 

u 2 = u 1 -. (91) 

K — 1 

Here, k is the polytropic index. For an infinitely thin shock wave front the divergence of the 
velocity is equal to 

div u = (wi -u 2 )S(X). (92) 

Substituting this expression into Eq. (87), we obtain that the charged particle acceleration 
at the fronts of collisionless shock waves propagating in a turbulent plasma is described by 
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V 



FIG. 19: Particle diffusion at the front of a shock wave propagating in a turbulent plasma. Inset: 
the phase plane of the particle bouncing between two plates. 

the equation (see Ref. [14] and references therein) 



In the limit, when the energy losses are negligibly small, this equation has a solution, which 
gives a power law dependence of the distribution function, / oc p~ k with the index value 
k = 3m2/(«2 — For k = 5/3 the index equals k = 4, i.e. / cx p~ A , or the energy spectrum 
dNcn{£) / d£ oc is close to the power law index observed in the galactic cosmic ray 
energy spectrum (see Eq. (83)). 

For the cosmic ray electron component in the high energy limit, at the energy when we 
cannot neglect the Compton and synchrotron losses, there is a cut off in the spectrum [153]. 
For typical parameters in supernova remnants, D = 10 25 cm 2 s _1 , B = 10 _4 G, u\ = 10 8 cm 
s^ 1 , the radiation losses limit the energy of ultrarelativistic electrons by values of the order 



Under the conditions of typical timescale of the laser plasmas the synchrotron losses of 
ultrarelativistic electrons interacting with the self-generated magnetic field is of the order of 



d x (u(X)f - Dd x f) + \d p (p 2 K(p)f) 



-2—^-5(X)\d p (p*f). 
3{k + 1) p l 



(93) 



of 10 TeV. 




(94) 
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IX. CONCLUSIONS 



Finally, we note that the development of superintense lasers with parameters in the ELI 
range will provide the necessary conditions for experimental physics where it will become 
possible to study ultrarelativistic energy of accelerated charged particles, super high intensity 
EMW and the relativistic plasma dynamics. A fundamental property of the plasma to create 
nonlinear coherent structures, such as relativistic solitons and vortices, collisionless shock 
waves and high energy particle beams, and to provide the conditions for relativistic regimes 
of the magnetic field line reconnection, makes the area of relativistic laser plasmas attractive 
for modeling of processes of key importance for relativistic astrophysics. 
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